Interactions of renormalized waves in thermalized Fermi-Pasta-Ulam chains 
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The dispersive interacting waves in Fermi-Pasta-Ulam (FPU) chains of particles in thermal equi- 
librium are studied from both statistical and wave resonance perspectives. It is shown that, even 
in a strongly nonlinear regime, the chain in thermal equilibrium can be effectively described by 
a system of weakly interacting renormalized nonlinear waves that possess (i) the Rayleigh- Jeans 
distribution and (ii) zero correlations between waves, just as noninteracting free waves would. This 
renormalization is achieved through a set of canonical transformations. The renormalized linear 
dispersion of these renormalized waves is obtained and shown to be in excellent agreement with nu- 
merical experiments. Moreover, a dynamical interpretation of the renormalization of the dispersion 
relation is provided via a self-consistency, mean-field argument. It turns out that this renormaliza- 
tion arises mainly from the trivial resonant wave interactions, i.e., interactions with no momentum 
exchange. Furthermore, using a multiple time-scale, statistical averaging method, we show that 
the interactions of near-resonant waves give rise to the broadening of the resonance peaks in the 
frequency spectrum of renormalized modes. The theoretical prediction for the resonance width for 
the thermalized /3-FPU chain is found to be in very good agreement with its numerically measured 
value. 



I. INTRODUCTION 

The study of discrete one-dimensional chains of par- 
ticles with the nearest-neighbor interactions provides in- 
sight to the dynamics of various physical and biological 
systems, such as crystals, wave systems, and biopoly- 
mers [3, [H, In the thermal equilibrium state, such 
nonlinear chains can be described by the canonical Gibbs 
measure U$ with the Hamiltonian 



(q- 3 -Q-j+i) 2 



+ V(qj -Qj+i), 



(1) 



where pj and qj are the momentum and the displace- 
ment from the equilibrium position of the j-th particle, 
respectively, V(qj — qj+i) is the anharmonic part of the 
potential, and the mass of each particle and the linear 
spring constant are scaled to unity. In this article, we 
only consider the potentials of the restoring type, i.e., the 
potentials for which the Gibbs measure exists. In order 
to study interactions of waves in such systems, one usu- 
ally introduces the canonical complex normal variables 
a k via 



a-k 



P k - lUJkQk 



(2) 



where Pk and Qk are the Fourier transforms of pj and qj, 
respectively, and u>k = 2 sin(7rfc/./V) is the linear disper- 
sion relation of the waves represented by a k . In terms of 
the ak, the Hamiltonian {1} becomes 



H 



E 



Wfc \ak\ 



V(a), 



(3) 



where V{a) is the combination of various products of a k 
and at corresponding to various wave-wave interactions. 



If the potential in Eq. {1} is harmonic, i.e., V = 0, then 
ak correspond to ideal, free waves, which have no energy 
exchanges among different k modes. In thermal equilib- 
rium, the Boltzmann distribution exp(— 0~ l ~^2ujk\a k \ 2 ) 
with temperature 0, gives rise to the following properties 
of free waves 



(a* k a t ) = n k 5f, 
= 0, 



[akai) 



(4) 
(5) 

for any k and I, where Uk = (|flfc| 2 ) = is the power 
spectrum. If the anharmonic part of the potential is 
sufficiently weak, then corresponding waves ak remain 
almost free, and Eqs. ((H) and §5§ would be approxi- 
mately satisfied in the weakly nonlinear regime. How- 
ever, when the nonlinearity becomes stronger, waves a k 
become strongly correlated, and, in general, the correla- 
tions between waves [Eq. ((SJ)] no longer vanish. In partic- 
ular, {a k aN-k) 7^ 0, as will be shown below. Naturally, 
the question arises: can the strongly nonlinear system 
in thermal equilibrium still be viewed as a system of al- 
most free waves in some statistical sense? In this article, 
we address this question with an affirmative answer: it 
turns out that the system §Q can be described by a com- 
plete set of renormalized canonical variables a k , which 
still possess the wave properties given by Eqs. ^ and (O 
with a renormalized linear dispersion. The waves that 
correspond to these new variables dk will be referred to 
as renormalized waves. Since these renormalized waves 
possess the equilibrium Rayleigh- Jeans distribution Q 
and vanishing correlations between waves, they resemble 
free, non-interacting waves, and can be viewed as statisti- 
cal normal modes. Furthermore, it will be demonstrated 
that the renormalized linear dispersion for these renor- 
malized waves has the form Cj k — r)(k)u>k, where r](k) is 
the linear frequency renormalization factor, and is inde- 
pendent of A: as a consequence of the Gibbs measure. 
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In our method, the construction of the renormalized 
variables dk does not depend on a particular form or 
strength of the anharmonic potential, as long as it is of 
the restoring type with only the nearest neighbor inter- 
actions, as in Eq. (fTJ). Therefore, our approach is non- 
perturbative and can be applied to a large class of sys- 
tems with strong nonlinearity. However, in this article, 
we will focus on the /3-FPU chain to illustrate the the- 
oretical framework of the renormalized waves. We will 
verify that dk effectively constitute normal modes for the 
/3-FPU chain in thermal equilibrium by showing that (i) 
the theoretically obtained renormalized linear dispersion 
relationship is in excellent agreement with its dynamical 
manifestation in our numerical simulation, and (ii) the 
equilibrium distribution of dk is still a Rayleigh-Jeans 
distribution and a^'s are uncorrelated. Note that similar 
expressions for the renormalization factor ry have been 
previously discussed in the framework of an approximate 
virial theorem [6| or effective long wave dynamics via the 
Zwanzig-Mori projection However, in our theory, the 
exact formula for the renormalization factor is derived 
from a precise mathematical construction of statistical 
normal modes, and is valid for all wave modes k — no 
longer restricted to long waves. 

Next, we address how renormalization arises from the 
dynamical wave interaction in the /3-FPU chain. We will 
show that the /3-FPU chain can be effectively described as 
a four- wave interacting Hamiltonian system of the renor- 
malized resonant waves dk- We will study the resonance 
structure of the /3-FPU chain and find that most of the 
exact resonant interactions are trivial, i.e., the interac- 
tions with no momentum exchange among different wave 
modes. In what follows, the renormalization of the linear 
dispersion will be explained as a collective effect of these 
trivial resonant interactions of the renormalized waves 
dk ■ We will use a self-consistency argument to find an ap- 
proximation, rise, of the renormalization factor r\. As will 
be seen below, the self-consistency argument essentially 
is of a mean-field type, i.e., the renormalization arises 
from the scattering of a wave by a mean-background of 
waves in thermal equilibrium via trivial resonant interac- 
tions. We note that our self-consistency, mean-field argu- 
ment is not limited to the weak nonlinearity. Very good 
agreement of the renormalization factor r\ and its dynam- 
ical approximation r\ sc — for weakly as well as strongly 
nonlinear waves — confirms that the renormalization is, 
indeed, a direct consequence of the trivial resonances. 

We will further study the properties of these renormal- 
ized waves by investigating how long these waves are co- 
herent, i.e., what their frequency widths are. Therefore, 
we consider near-resonant interactions of the renormal- 
ized waves dk, i.e., interactions that occur in the vicinity 
of the resonance manifold, since most of the exact res- 
onant interactions are trivial, i.e., with no momentum 
exchanges, and they, cannot effectively redistribute en- 
ergy among the wave modes. 

We will demonstrate that near-resonant interactions 
of the renormalized waves dk provide a mechanism for 



effective energy exchanges among different wave modes. 
Taking into account the near-resonant interactions, we 
will study analytically the frequency peak broadening of 
the renormalized waves dk by employing a multiple time- 
scale, statistical averaging method. Here, we will arrive 
at a theoretical prediction of the spatiotemporal spec- 
trum |afc(aj)| 2 , where afe(w) is the Fourier transform of 
the normal variable 5fc(t), and u> is the frequency. The 
predicted width of frequency peaks is found to be in good 
agreement with its numerically measured values. 

In addition, for a finite /3-FPU chain, we will mention 
the consequence, to the correlation times of waves, of the 
momentum exchanges that cross over the first Brillouin 
zone. This process is known as the umklapp scattering 
in the setting of phonon scattering Note that, in 
the previous studies @ of the FPU chain from the wave 
turbulence point of view, the effects arising from the finite 
nature of the chain were not taken into account, i.e., only 
the limiting case of N — > oo, where N is the system size, 
was considered. 

The article is organized as follows. In Section [H] we 
discuss a chain of particles with the nearest-neighbor non- 
linear interactions. We demonstrate how to describe a 
strongly nonlinear system as a system of waves that re- 
semble free waves in terms of the power spectrum and 
vanishing correlations between waves. We show how to 
construct the corresponding renormalized variables with 
the renormalized linear dispersion. In Section IIII1 we 
rewrite the /3-FPU chain as an interacting four-wave 
Hamiltonian system. We study the dynamics of the 
chain numerically and find excellent agreement between 
the renormalized dispersion, obtained analytically and 
numerically. In Section HVi we describe the resonance 
manifold analytically and illustrate its controlling role 
in long-time averaged dynamics using numerical simula- 
tion. In Section [Vj we derive an approximation for the 
renormalization factor for the linear dispersion using a 
self-consistency condition. In Section IVH we study the 
broadening effect of frequency peaks and predict analyt- 
ically the form of the spatiotemporal spectrum for the 
/3-FPU chain. We provide the comparison of our pre- 
diction with the numerical experiment. We present the 
conclusions in Section IVII1 



II. RENORMALIZED WAVES 

Consider a chain of particles coupled via nonlinear 
springs. Suppose the total number of particles is N and 
the momentum and displacement from the equilibrium 
position of the j-th particle are pj and qj, respectively. 
If only the nearest-neighbor interactions are present, then 
the chain can be described by the Hamiltonian 

H = H 2 + V, (6) 
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where the quadratic part of the Hamiltonian takes the 
form 



Hi 



N 
\ " 2 



3 = 1 



(7) 



and the anharmonic potential V is the function of the 
relative displacement qj — qj+%. Here periodic boundary 
conditions qn+i = 9i and pn+i = Pi are imposed. Since 
the total momentum of the system is conserved, it can 
be set to zero. 

In order to study the distribution of energy among the 
wave modes, we transform the Hamiltonian to the Fourier 
variables Qk, Pk via 



N-l 



27rifcj 

q je w , 

3=0 
1 N ~ 1 



m 



This transformation is canonical 
tonian ([6|) becomes 



10, 111 and the Hamil- 



H 



JV-l 

2 ^ 

k=l 



P k \ 2 +cj 2 \Q k \ 2 + V(Q), 



(9) 



If the nonlinear interactions are weak, then it is con- 
venient to further transform the Hamiltonian © to the 
complex normal variables defined by Eq. This trans- 
formation is canonical, i.e., the dynamical equation of 
motion becomes 



ia k 



dH 



(13) 



In terms of these normal variables, the Hamiltonian © 
takes the form ([3]). For the system of noninteracting 

waves, i.e., H = 53fc=i w fcl a k| 2 ; we obtain a standard 
virial theorem in the form 



(K k )\ v=0 = (U k )\ v =o- 



(14) 



As a consequence of this virial theorem, we have the prop- 
erties of free waves, which were already mentioned above 
[Eqs. and ©], i.e., 

(a%ai) = ^-{{\P k \ 2 )+u3 2 {\Q k \ 2 ))5^-8l (15) 

ZU) k LU k 

M = ^((\P k \ 2 )-u; 2 (\Q k \ 2 ))6 k N +l =0, (16) 



for all wave numbers k and I. Note that equation (fl5() 
gives the classical Rayleigh- Jeans distribution for the 
power spectrum of free waves 



where to k = 2sm(nk/N) is the linear dispersion relation. 
Note that, throughout the paper, for the simplicity of 
notation, we denote the periodic wave number space by 
the set of integers in the range [0, N — 1], i.e., we drop 
the conventional factor, 2ir/N. The zeroth mode vanishes 
due to the fact that the total momentum is zero. 

If the system © is in thermal equilibrium, then the 
canonical Gibbs measure, with the corresponding parti- 
tion function 



(10) 



with the temperature 9, can be used to describe the sta- 
tistical behavior of the system. We consider the systems 
with the anharmonic potential of the restoring type, i.e., 
the potential for which the integral in Eq. (|10p converges. 
It can be easily shown that for system ([9]) the average ki- 
netic energy K k of each mode is independent of the wave 
number 



{Kk) = {Ki), 



(11) 



where k and I are wave numbers, K k = \P k \ 2 /2, and (. . .) 
denotes averaging over the Gibbs measure. Similarly, the 
average quadratic potential U k of each mode is indepen- 
dent of the wave number 



(Uk) = (Ui), 



(12) 



where U k = tu 2 \Q k \ 2 /2. 



n k = — . 



(17) 



However, if the nonlinearity is present, the waves a k and 
ajv_fc become correlated, i.e., 



(a k a N - k ) 



1 

2uj, 



((\P k \ 2 )-u 2 (\Q k \ 2 ))^0, (18) 



since the property (jT4j) is no longer valid. 

As we mentioned before, a complete set of new renor- 
malized variables a k can be constructed, so that the 
strongly nonlinear system can be viewed as a system of 
"free" waves in the sense of vanishing correlations and 
the power spectrum, i.e., the new variables a k satisfy 
the properties of free waves given in Eqs. (TTS")) and (|16p . 
Next, we show how to construct these renormalized vari- 
ables a k . 

Consider the generalization of the transformation ((2|) , 
namely, the transformation from the Fourier variables Q k 
and P k to the renormalized variables a k by 



(ik 



P k - \u} k Q k 



(19) 



where u>k is an arbitrary function with the only restric- 
tions 



(20) 



One can show that, these restrictions (f2"D|) provide a nec- 
essary and sufficient condition for the transformation (TT 
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to be canonical. For the renormalized waves a,k, we can 
compute 

<a*5i> = i^- k {(\Pk\ 2 )+ul{\Q k \ 2 ))5l (21) 

(5fc5,) = ^ k ((\ P ^)-^t(\Q^)) 6 N +l - ( 22 ) 

Since we have the freedom of choosing any uj k (with 
the only restrictions (|2"0|) ). we can chose u> k such that 
(dkCLN-k) vanishes. Thus, the renormalized variables a k 
for a strongly nonlinear system will behave like the bare 
variables a% for a noninteracting system in terms of van- 
ishing correlations between waves. Therefore, we deter- 
mine u) k via 



(l^l 2 ) 



-u 2 k (\Q k \ 2 )=Q. 



(23) 



Note that the requirement (|23[) has the form of the virial 
theorem for the free waves but with the renormalized 
linear dispersion ui k . We rewrite Eq. (|2"3")) in terms of the 
kinetic and quadratic potential parts of the energy of the 
mode k as 



(24) 



The in Eqs. (fTTjl and (| 1 2[) leads to the k independence of 
the right-hand side of Eq. (|24p . This allows us to define 
the renormalization factor r\ for all k's by 




= u^_ (K)_ 
V ~ V (u) ' 



(25) 



for dispersion tOk- Here K 



^2 k=1 Uk are the kinetic and the quadratic potential 
parts of the total energy of the system © , respectively. 
Note that the way of constructing the renormalized vari- 
ables Ofc via the precise requirement of vanishing corre- 
lations between waves yields the exact expression for the 
renormalization factor, which is valid for all wave num- 
bers k and any strength of nonlinearity. The indepen- 
dence of rj of the wave number k is a consequence of the 
Gibbs measure. This k independence phenomenon has 
been observed in previous numerical experiments [U, [l2T ] . 
We will elaborate on this point in the results of the nu- 
merical experiment presented in Section IIIII 

The immediate consequence of the fact that r\ is inde- 
pendent of k is that the power spectrum of the renormal- 
ized waves possesses the precise Rayleigh- Jeans distribu- 
tion, i.e., 



n k = — , 



(26) 



from Eq. (|2"Tj) . where h k = (|afc| 2 ). Combining Eqs. ([2]) 
and (fT9|) . we find the relation between the "bare" waves 
a k and the renormalized waves a k to be 

ak = \(^ + ^)~ ak + \(^ ] --^~ aN - k - (2?) 



Using Eq. (|27p . we obtain the following form of the power 
spectrum for the bare waves a k 



If 1 \ 9 

2 V T J Uk 



(28) 



which is a modified Rayleigh- Jeans distribution due to 
the renormalization factor (1 + l/rj 2 )/2. Naturally, if the 
nonlinearity becomes weak, we have rj — ► 1, and, there- 
fore, all the variables and parameters with tildes reduce 
to the corresponding "bare" quantities, in particular, 
Lu k — > uj k , a k — > a k , h k — > n k . It is interesting to point 
out that, even in a strongly nonlinear regime, the "free- 
wave" form of the Rayleigh- Jeans distribution is satisfied 
exactly [Eq. (|26[) ] by the renormalized waves. Thus, we 
have demonstrated that even in the presence of strong 
nonlinearity, the system in thermal equilibrium can still 
be viewed statistically as a system of "free" waves in the 
sense of vanishing correlations between waves and the 
power spectrum. 

Note that, in the derivation of the formula for the 
renormalization factor [Eq. (|2"5j) ]. we only assumed the 
nearest-neighbor interactions, i.e., the potential is the 
function of qj — Qj+i- One of the well-known examples of 
such a system is the /3-FPU chain, where only the forth 
order nonlinear term in V is present. In the remainder of 
the article, we will focus on the /3-FPU to illustrate the 
framework of the renormalized waves a k . 



III. NUMERICAL STUDY OF THE /3-FPU 
CHAIN 



Since its introduction in the early 1950s, the study of 
the FPU lattice [13| has led to many great discoveries 
in mathematics and physics, such as soliton theory [3[. 
Being non-integrable, the FPU system also became in- 
tertwined with the celebrated Kolmogorov-Arnold-Moser 
theorem [11| . Here, we extend our results of the thermal- 
ized /3-FPU chain, which were briefly reported in (l2| . 

The Hamiltonian of the /3-FPU chain is of the form 



N 1 1 /3 



(29) 



3=1 



where /3 is a parameter that characterizes the strength of 
nonlinearity. 

The canonical equations of motion of the /3-FPU chain 
are 



dH 
dH 



(30) 



^ = ~ dq~ = ~ 2<li + qj+1 ^ 

+ 0[(.Qj+i ~ ft') 3 - (Qj ~ Ij-lf]- 

To investigate the dynamical manifestation of the renor- 
malized dispersion Cu k of a k , we numerically integrate 
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Eq. pOl) . Since we study the thermal equilibrium 
state [3 EE EE E3] of the /3-FPU chain, we use ran- 
dom initial conditions, i.e., pj and qj are selected at 
random from the uniform distribution in the intervals 
(-PmaxjPmax) and (- g max , g max ) , respectively, with the 
two constraints that (i) the total momentum of the sys- 
tem is zero and (ii) the total energy of the system E 
is set to be a specified constant. We have verified that 
the results discussed in the paper do not depend on de- 
tails of the initial data. Note that the behavior of (3- 
FPU for fixed TV is fully characterized by only one pa- 
rameter (3E [l8|. We use the sixth order symplectic 
Yoshida algorithm with the time step dt = 0.01, 
which ensures the conservation of the total system en- 
ergy up to the ninth significant digit for a runtime 
t = 10 6 time units. In order to confirm that the sys- 
tem has reached the thermal equilibrium state [2p| . the 
value of the energy localization [2lj was monitored via 

L(t) = JVXjLi Gy/CEjLi Gj) 2 , where G s is the energy 
of the j-th particle defined as 



G j = \v) + \ [(«," - Qj+i) 2 + {Qj- 



i - Qjj 

41 



(31) 



If the energy of the system is concentrated around one 
site, then L(t) — O(N). Whereas, if the energy is uni- 
formly distributed along the chain, then L(t) = 0(1). 
In our simulations, in thermal equilibrium states, L(t) is 
fluctuating in the range of 1-3. Since our simulation is 
of microcanonical ensemble, we have monitored various 
statistics of the system to verify that the thermal equi- 
librium state that is consistent with the Gibbs distribu- 
tion (canonical ensemble) has been reached. Moreover, 
we verified that, for TV as small as 32 and up to as large 
as 1024, the equilibrium distribution in the thermalized 
state in our microcanonical ensemble simulation is con- 
sistent with the Gibbs measure. We compared the renor- 
malization factor ([25jl by computing the values of (K) 
and (U) numerically and theoretically using the Gibbs 
measure and found the discrepancy of r\ to be within 
0.1% for {3 = 1 and the energy density E/N — 0.5 for TV 
from 32 to 1024. 

We now address numerically how the renormalized lin- 
ear dispersion uik manifests itself in the dynamics of the 
/3-FPU system. We compute the spatiotemporal spec- 
trum |afc(tj)| 2 , where &k(oj) is the Fourier transform of 
5fc(t). (Note that, for simplicity of notation, we drop a 
tilde in dk) Figure [1] displays the spatiotemporal spec- 
trum of ctfc, obtained from the simulation of the /3-FPU 
chain for TV = 256, (3 = 0.5, and E = 100. In order to 
measure the value of r\ from the spatiotemporal spectrum, 
we use the following procedure. For the fixed wave num- 
ber k, the corresponding renormalization factor r)(k) is 
determined by the location of the center of the frequency 
spectrum |a / t(cj)|' ; 



i.e. 



T}(k) = 



^c(fc) 



with u> c (k) = 



fuj\a k (uj)\ 2 duo 
J \a k (uj)\ 2 doj 




0.5 1 



FIG. 1: The spatiotemporal spectrum |afc(aj)| 2 in thermal 
equilibrium. The chain was modeled for TV = 256, (3 — 0.5, 
and E = 100. [max{— 8, In |ofe(tj)| 2 }, with corresponding gray 
scale, is plotted for a clear presentation]. The solid curve 
corresponds to the usual linear dispersion u)k = 2sin(7rfc/TV). 
The dashed curve shows the locations of the actual frequency 
peaks of \ak(io)\ 2 . 



The renormalization factor rj{k) of each wave mode k is 
shown in Fig. [2] (inset). The numerical approximation fj 
to the value of rj is obtained by averaging all r)(k), i.e., 



1 



N-l 



N- 



fc=i 



The renormalization factor for the case shown in Fig. [T] 
is measured to be fj ss 1.1824. It can be clearly seen in 
Fig. [2] (inset) that ry(fc) is nearly independent of k and 
its variations around fj are less than 0.3%. We also com- 
pare the renormalization factor rj obtained from Eq. (|25|) 
(solid line in Fig. [2] (inset)) with its numerically com- 
puted approximation fj (dashed line in Fig. [2] (inset)). 
Equation ([23)) gives the value rj rj 1.1812 and the differ- 
ence between r\ and fj is less then 0.1%, which can be 
attributed to the statistical errors in the numerical mea- 
surement. In Fig. [51 we plot the value of rj as a function 
of (3 for the system with TV = 256 particles and the to- 
tal energy E = 100. The solid curve was obtained using 
Eq. (|25p while the circles correspond to the value of T) 
determined via the numerical spectrum |afc(ix>)| 2 as dis- 
cussed above. It can be observed that there is excellent 
agreement between the theoretic prediction and numeri- 
cally measured values for a wide range of the nonlinearity 
strength (3. 

In the following Sections, we will discuss how the renor- 
malization of the linear dispersion of the /3-FPU chain in 
thermal equilibrium can be explained from the wave res- 
onance point of view. In order to give a wave description 
of the /3-FPU chain, we rewrite the Hamiltonian ([2U)l in 
terms of the renormalized variables [Eq. (|19p ] with 
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FIG. 2: The renormalization factor as a function of the non- 
linearity strength f3. The analytical prediction [Eq. (|25[) ] is 
depicted with a solid line and the numerical measurement is 
shown with circles. The chain was modeled for N = 256, and 
E = 100. Inset: Independence of k of the renormalization 
factor n(k). The circles correspond to n(k) obtained from the 
spatiotemporal spectrum shown in Fig. [1] [only even values 
of k are shown for clarity of presentation]. The dashed line 
corresponds to the mean value fj. For /3 = 0.5, the mean value 
of the renormalization factor is found to be fj ~ 1.1824. The 
variations of rjk around fj are less then 0.3%. [Note the scale of 
the ordinate.] The solid line corresponds to the renormaliza- 
tion factor n obtained from Eq. (|25[1 ■ For the given parameters 
n » 1.1812. 



Wfe = TjUJ k , 



N-l , . 



a/c| 2 



fc=i 

+ "J" \jl - ^ ) I '*'/,'*' A -A ' 



N-l 



2 



+ [ -Af m a k &i& m a*+c.c. 



+ ( -A$ lms a k aia m a s + c.c. 



where c.c. stands for complex conjugate, and 



(32) 



-A7 



3/3 
SNrJ- 



■y/LUkUJlUl m UJs 



(33) 



is the interaction tensor coefficient. Note that, due to the 
discrete nature of the system of finite size, the wave space 
is periodic and, therefore, the "momentum" conservation 
is guaranteed by the following "periodic" Kronecker delta 
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(34) 
(35) 
(36) 



Here, the Kronecker (^-function is equal to 1, if the sum 
of all superscripts is equal to the sum of all subscripts, 
and 0, otherwise. 



IV. DISPERSION RELATION AND 
RESONANCES 

In order to address how the renormalized dispersion 
arises from wave interactions, we study the resonance 
structure of our nonlinear waves. Since the system (|32|) 
is a Hamiltonian system with four-wave interactions, we 
will discuss the properties of the resonance manifold asso- 
ciated with the /3-FPU system described by Eq. (|32| as 
a first step towards the understanding of its long time 
statistical behavior. We comment that the resonance 
structure is one of the main objects of investigation in 
wave turbulence theory @, [H [23|, H, [H [H, E3] ■ The 
theory of wave turbulence focuses on the specific type of 
interactions, namely resonant interactions, which domi- 
nate long time statistical properties of the system. On 
the other hand, the non-resonant interactions are usually 
shown to have a total vanishing average contribution to 
a long time dynamics. 

In analogy with quantum mechanics, where a + and a 
are creation and annihilation operators, we can view a£ 
as the outgoing wave with frequency Cbk and dk as the 
incoming wave with frequency ibk- Then, the nonlinear 
term a£a*a m ci s A^ s in system (f3"2j) can be interpreted as 
the interaction process of the type (2 — > 2), namely, two 
outgoing waves with wave numbers k and / are "created" 
as a result of interaction of the two incoming waves with 
wave numbers m and s. Similarly, a^aia m a s A i ^ !ls in sys- 
tem (l3"2"j) describes the interaction process of the type 
(3 — * 1), that is, one outgoing wave with wave number k 
is "created" as a result of interaction of the three incom- 
ing waves with wave numbers I, m, and s, respectively. 
Finally, dkaia m a s AQ lms describes the interaction process 
of the type (4 — > 0), i.e., all four incoming waves inter- 
act and annihilate themselves. Furthermore, the complex 



conjugate terms 5^5,,* a^a 



■ A k 

■ ^hns 



and a* k d*a^a 



"A Hp- 



3) and 



scribe the interaction processes of the type (1 
(0 — ► 4), respectively. 

Instead of the processes with the "momentum" conser- 
vation given via the usual <5^ s , 5^ lm , or Sq 1 " 111 functions 
for an infinite discrete system, the resonant processes of 
the /3-FPU chain of a finite size are constrained to the 



manifold given by A„' s , Af 171 , or 



A klms 
^0 ! 



respectively. 



Next, we describe these resonant manifolds in detail. As 
will be pointed out in Section IVTl there is a consequence 
of this finite size effect to the properties of the renormal- 
ized waves. 
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The resonance manifold that corresponds to the (2 — > 
2) resonant processes in the discrete periodic system, 
therefore, is described by 



N 



k + 1 = m + s, 
ij k + Cji — Q m + Q s , 



(37) 



/V 



where we have introduced the notation g = h, which 
means that g — h, g = h + N , or g — h — N for any g 
and h. The first equation in system (|37[) is the "momen- 
tum" conservation condition in the periodic wave num- 
ber space. This "momentum" conservation comes from 
I J = 1- (Note that |A*J S | can assume only the value 
of 1 or 0.) Similarly, from \Af m \ = 1 and |A(f ms | = 1, 
the resonance manifolds corresponding to the resonant 
processes of types (3 — > 1) and (4 — * 0) are given by 



and 



{k + I + to = s, 
Wfc + a)/ + £> m = a) s , 

k + l + m + s = 0, 

U>k + + &m + w s = 0, 



(38) 



(39) 



respectively. For the processes of type (3 — » 1), the nota- 
tion g = h means that g — h, g = h + N, or g^h + 2N . 

For the (4 — > 0) processes, g = h means that g = h + N , 
.g = h + 2N, or .g = h + 3N. 

To solve system (|37j). we rewrite it in a continuous 
form with x = fe/iV, y = l/N , z — m/N, v — s/N , which 
are real numbers in the interval (0, 1). By recalling that 
Ciik = 2t) sin(n k/N) , we have 



x + y = z + v, 

sin(7rx) + sin(7ry) = sin(7rz) + sin(7ru) 



(40) 



Thus, any rational quartet that satisfies Eq. (I40p yields 
a solution for Eq. (|37|) . There are two distinct types of 
the solutions of Eq. J30|). The first one corresponds to 
the case 

x + y = z + v, 
whose only solution is given by 



(41) 



i.e., these are trivial resonances, as we mentioned above. 
The second type of the resonance manifold of the (2 — > 
2)-type interaction processes corresponds to 

x + y = z + v±l, 

the solution of which can be described by the following 
two branches 





Z-2 



x + y 1 

2 7T 

x + y 



arcsin(A) + 2j, 



(42) 



1 arcsin(A) + 2j, (43) 

7T 



where A = tan (tt(x + y)/2) cos (tt(x — y)/2) and j is an 
integer. The second type of resonances arises from the 
discreteness of our model of a finite length, leading to 
non-trivial resonances. For our linear dispersion here, 
non-trivial resonances are only those resonances that in- 
volve wave numbers crossing the first Brillouin zone. As 
mentioned above, in the setting of the phonon physics, 
these non-trivial resonant processes are also known as 
the umklapp scattering processes. In Fig. [3l we plot the 
solution of Eq. (|40|) for x = k/N with the wave num- 
ber k — 90 for the system with N = 256 particles (the 
values of k and N are chosen merely for the purpose 
of illustration). We stress that all the solutions of the 
system ((40j) are given by the Eqs. ijil]). (l42|) . and (|43|). 
and that the non-trivial solutions arise only as a conse- 
quence of discreteness of the finite chain. The curves in 
Fig. [3] represent the loci of (z,y), parametrized by the 
fourth wave number v, i.e., x, y, z, and v form a reso- 
nant quartet, where z — m/N, and y — l/N. Note that 
the fourth wave number v is specified by the "momen- 
tum" conservation, i.e., the first equation in Eq. (I40p . 
The two straight lines in Fig. [3] correspond to the trivial 
solutions, as given by Eq. (|4"Tj) . The two curves (dot- 
ted and dashed) depict the non-trivial resonances. Note 
that the dotted part of non-trivial resonance curves cor- 
responds to the branch (j42|) . and the dashed part corre- 
sponds to the branch (|4"3"]) , respectively. An immediate 
question arises: how do these resonant structures man- 
ifest themselves in the FPU dynamics in the thermal 
equilibrium? By examining the Hamiltonian (|32|) . we 
notice that the resonance will control the contribution of 
terms like a^a*a m a s Aj^ s in the long time limit. There- 
fore, we address the effect of resonance by computing 
long time average, i.e., (a£a*a m a s )A^,., and comparing 
this average (Fig. gj with Fig. El To obtain Fig. H the 
/3-FPU system was simulated with the following param- 
eters: N — 256, (3 — 0.5, E — 100, and the averaging 
time window r = 400ti, where t\ is the longest linear 
period, i.e., t\ = 2-k/H)\. In Fig. [U mode k was fixed 
with k = 90 and the mode s, a function of k 7 I, and 

m, is obtained from the constraint k + I = m + s, i.e., 
|Am S | = 1. Note that we do not impose here the condi- 
tion (Dfc = G) m +Q S , therefore, \ (a^a* a m a s ) A^[ s \ is a 
function of I and m. By comparing Figs.[3]and|31 it can be 
observed that the locations of the peaks of the long time 
average |(a£a ; *o m a s )Aj^ s | coincide with the loci of the 
(2 — > 2)-type resonances. This observation demonstrates 
that, indeed, there are nontrivial (2 — > 2)-type resonances 
in the finite /3-FPU chain in thermal equilibrium. Fur- 
thermore, it can be observed in Fig.g]that, in addition to 
the fact that the resonances manifest themselves as the 
locations of the peaks of \{a^.a*a m a s )A^ s \, the structure 
of near-resonances is reflected in the finite width of the 
peaks around the loci of the exact resonances. Note that, 
due to the discrete nature of the finite /3-FPU system, 
only those solutions x, y, z, and v of Eq. (|40[) . for which 
Nx, Ny, Nz, and Nv are integers, yield solutions k, I, 
to, and s for Eq. ([3T]) . In general, the rigorous treatment 
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z 



FIG. 3: The solutions of Eq. p)]). The solid straight lines 
correspond to the trivial resonances [solutions of Eq. (|41[) ]. 
The solutions are shown for fixed x = k/N, k — 90, N — 256 
as the fourth wave number v scans from 1/7V to (N — 1)/N 
in the resonant quartet Eq. (|40p . The non-trivial resonances 
are described by the dotted or dashed curves. The dotted 
branch of the curves corresponds to the non-trivial resonances 
described by Eq. (I42p and and the dashed branch corresponds 
to the non-trivial resonances described by Eq. (|43p . 




0.2 0.4 0.6 0.8 1 

m/N 



FIG. 4: The long time average j (a* k a*a m a s ) A^ s | of the /3-FPU 
system in thermal equilibrium. The parameters for the FPU 
chain are N = 256, (3 = 0.5, and E = 100. (ala*dma s ) A^ s 
was computed for fixed k = 90. The darker grayscale cor- 
responds to the larger value of (a^a*a rn a s )A^ 1 [ 8 . The exact 
solutions of Eq. (|40p . which are shown in Fig. [3] coincide 
with the locations of the peaks of |{5£S^S m S s )A(^ s |. There- 
fore, the darker areas represent the near-resonance struc- 
ture of the finite /3-FPU chain. (The two white lines show 
the locations, where s = and, therefore, a£5 ; *a m 5 s Aj^ s = 
0.) [max{2, ln(j (a£a;*5 m 5 s ) Aj^ s |)} with the corresponding 
grayscale is plotted for a clean presentation]. 



of the exact integer solutions of Eq. (f37|) is not straight- 
forward. For example, for N — 256, we have the follow- 
ing two exact quartets k — {k, l,m, s}: k = {k,N/2 — 
k, N/2 + k,N - k}, k = {k, N/2 -k,N - k, N/2 + k} 
for k < N/2, and k = {k, 3N/2 - k,k- N/2, N - k}, 
k = {k,3N/2~k,N-k,k-N/2} for k > N/2. We have 
verified numerically that for N = 256 there are no other 
exact integer solutions of Eq. (|37|) . In the analysis of the 
resonance width in Section IVT1 we will use the fact that 
the number of exact non-trivial resonances [Eq. ([37)) ] is 
significantly smaller than the total number of modes. 

The broadening of the resonance peaks in Fig. 2] sug- 
gests that, to capture the near-resonances for character- 
izing long time statistical behavior of the /3-FPU system 
in thermal equilibrium, instead of Eq. (|37|) , one needs to 
consider the following effective system 



\u>k + &l — & m - ui s \ < Alu, 

where < Aw <C LUk for any k, and Alu characterizes the 
resonance width, which results from the near-resonace 
structure. Clearly, Auj is related to the broadening of 
the spectral peak of each wave d a (t) with a = k,l,m, 
or s in the quartet, and this broadening effect will be 
studied in detail in Section [VlJ Note that the structure 
of near-resonances is a common characteristic of many 
periodic discrete nonlinear wave systems [HI, [29|, |30[ ■ 

Further, it is easy to show that the dispersion rela- 
tion of the /3-FPU chain does not allow for the occur- 
rence of (3 — > l)-type resonances, i.e., there are no solu- 
tions for Eq. (f3"5)) . and, therefore, all the nonlinear terms 
a^amOsA^ are non-resonant and their long time av- 
erage {a* k aia m a s ) A^ ms vanishes. As for the resonances 
of type (4 — > 0), since the dispersion relation is non- 
negative, one can immediately conclude that the solution 
of the system (l39l) consists only of zero modes. Therefore, 
the processes of type (4 — > 0) are also non-resonant, giv- 
ing rise to (akaia m a s )A(f ms = 0. In this article, we will 
neglect the higher order effects of the near-resonances of 
the types (3 -> 1) and (4 -> 0). 

In the following sections, we will study the effects of 
the resonant terms of type (2 — + 2), namely, the linear 
dispersion renormalization and the broadening of the fre- 
quency peaks of dk(t). It turns out that, the former is 
related to the trivial resonance of type (2 — > 2) and the 
latter is related to the near-resonances, as will be seen 
below. 



V. SELF-CONSISTENCY APPROACH TO 
FREQUENCY RENORMALIZATION 

We now turn to the discussion of how the trivial reso- 
nances give rise to the dispersion renormalization. This 
question was examined in [l2j before. There, it was 
shown that the renormalization of the linear dispersion 
of the /3-FPU chain arises due to the collective effect of 
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the nonlinearity. In particular, the trivial resonant inter- 
actions of type (2 — > 2), i.e., the solutions of Eq. (|4"T) . 
enhance the linear dispersion (the renormalized disper- 
sion relation takes the form ujk = r/LUk with 77 > 1), and 
effectively weaken the nonlinear interactions. Here, we 
further address this issue and present a self-consistency 
argument to arrive at an approximation for the renor- 
malization factor 7/. As it was mentioned above, the 
contribution of the non-resonant terms have a vanishing 
long time effect to the statistical properties of the sys- 
tem, therefore, in our self-consistent approach, we ignore 
these non-resonant terms. By removing the non-resonant 
terms and using the canonical transformation 



P fc - \r) sc u k Q k 



(45) 



where r) sc is a factor to be determined, we arrive at a sim- 
plified effective Hamiltonian from Eq. (f3"2")l for the finite 
/3-FPU system 



N-l 



CJfe 



1 

Vs, 



\d k \ 2 



(46) 



fc=i 

N-l 

+ 2^ T m S A ms a k a l a m a s- 

fc,i,m,s— 1 

The "off-diagonal" quadratic terms dkdN-k from 
Eq. (1321) are not present in Eq. (j46|) . since dk are chosen 
so that {a k aN-k) = (see Section [TTJ) . The contribution 
of the trivial resonances in H c g is 



N-l 



2 |a fc | 2 , 



(47) 



k,l=l 



which can be "linearized" in the sense that averaging the 
coefficient in front of |afc| 2 in iJ| r gives rise to a quadratic 
form 



H. 



N-l 
tr _ \^ 



N-l 



k=l 



x |2 



1=1 




Note that the subscript 2 in if| r emphasizes the fact that 
H% T now can be viewed as a Hamiltonian for the free 
waves with the familiar effective linear dispersion Qk = 

4 Efc7 lT fe/(N 2 > SHI]- This linearization is essentially 
a mean-field approximation, since the long-time average 
of trivial resonances in Eq. (|47[) is approximated by the 
interaction of waves dk with background waves (|a;|). 
The self consistency condition, which determines rj sc , can 
be imposed as follows: the quadratic part of the Hamil- 
tonian (|46[) . combined with the "linearized" quadratic 
part, H2, of the quartic Hf, should be equal to an effec- 
tive quadratic Hamiltonian H2 = ^2k=i &k\a>k\ 2 for the 



renormalized waves, i.e., 



7V-1 



E 



Wk 



k=l 

N-l / N-l 



\~a k \ 2 



2 \ Vs 
£ (4^^/(h| 2 ))|a fc | 2 =^^|5 fc | 2 , 



N-l 



k=l \ 1 = 1 



k=l 



(48) 



where ibk is the renormalized linear dispersion, which is 
used in the definition of our renormalized wave, Eq. (|45[) . 
and ibk = r] sc tOk- Equating the coefficients of t^lofcl 2 on 
both sides for every wave number k yields 



\ 1 



N-l 



'Isc J ,_. 



3/3 



where use is made of Eq. (|3"3"|) . After algebraic simplifi- 
cation, we have the following equation for rj sc 



N-l 

'/:;, '/>.• : : — E ^(N 2 )- 
1=1 



3(3 

77 



(49) 



Using the property (|2ip of the renormalized normal vari- 
ables dk, we find the following dependence of (|flfc| 2 ) on 

Vsc, 



(N 



f 



(dp^+^fim 2 )). (50) 



Combining Eqs. (|4^|) and ([50]) leads to 

V 4 SC - A V 2 C - B = 0, (51) 

where 



;=i 



The only physically relevant solution of Eq. ([511 is 



1 A + y/A 2 + AB 



(52) 



The constants A and B can be easily derived using the 
Gibbs measure. 

Next, we compare the renormalization factor r\ 
[Eq. (|2"5|) ] with its approximation r/ sc [Eq. ([52"| ] from the 
self-consistency argument. In Appendix [XJ we study in 
detail the behavior of both r\ and r\ sc in the two limiting 
cases, i.e., when nonlinearity is small {[3 — > with fixed 
total energy E), and when nonlinearity is large {(3 — > 00 
with fixed total energy E). As is shown in Appendix YK[ 
for the case of small nonlinearity, both rj and rj sc have the 
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FIG. 5: The renormalization factor as a function of the non- 
linearity strength (3 for small values of /3. The renormalization 
factor r] [Eq. (|25[) ] is shown with the solid line. The approx- 
imation r] sc [Eq. (|52p ](via the self-consistency argument) is 
depicted with diamonds connected with the dashed line. The 
small-/? limit [Eq. (I53p ] is shown with the solid circles con- 
nected with the dotted line. Note that, abscissa is of loga- 
rithmic scale. Inset: The renormalization factor as a func- 
tion of the nonlinearity strength /3 for large values of /3. The 
renormalization factor r\ [Eq. (|25|l ] is shown with the solid 
line. r/ sc [Eq. (|52p ] is depicted with diamonds connected with 
the dashed line. The large-/? scaling [Eq. (|54|l ] is shown with 
the dashed-dotted line. Note that, the plot is of log-log scale 
with base 10. 



same asymptotic behavior in the first order of the small 
parameter /?, 



3B 

V = 1 + — /3 + 0(/? 2 ), 
Vsc = l + ^/3 + 0(/3 2 ). 



Moreover, in the case of strong nonlinearity (3 
r\ and rj sc scale as i.e., 

V ~ Vsc ~ (3^ 



(53) 



both 



(54) 



(see Appendix lAl for details). Note that, in [T2|, we nu- 
merically obtained the scaling 77 ~ P ' 2 , which differs 
from the exact analytical result (154|) due to statistical er- 
rors in the numerical estimate of the power. In Fig.0 we 
plot the renormalization factor 77 and its approximation 
ri sc for the case of small nonlinearity (3 for the system with 
N = 256 particles and total energy E = 100. The solid 
line shows 77 computed via Eq. (f2"5"]) . the diamonds with 
the dashed line represent the approximation via Eq. (|52p . 
and the solid circles with the dotted line correspond to 
the small-/? limit (fSU)) . In Fig. [5] (inset), we plot the renor- 
malization factor 77 and its approximation rj sc for the case 



of large nonlinearity (3 for the system with N = 256 par- 
ticles and total energy E = 100. The solid line shows 
77 computed via Eq. (|2"5"j). the diamonds with the dashed 
line represent the approximation via Eq. (|52p . and the 
dashed-dotted line correspond to the large-/3 scaling f|54[) . 
Figure O shows good agreement between the renormaliza- 
tion factor 77 and its approximation rj sc from the self con- 
sistency argument for a wide range of nonlinearity, from 
(3 ~ 10 -3 to (3 ~ 10 4 . This agreement demonstrates, that 
(i) the effect of the linear dispersion renormalization, in- 
deed, arises mainly from the trivial four-wave resonant 
interactions, and (ii) our self-consistency, mean- field ar- 
gument is not restricted to small nonlinearity. 



VI. RESONANCE WIDTH 

Finally, we address the question of how coherent these 
renormalized waves are, i.e., we study how the nonlin- 
ear interactions of waves in thermal equilibrium broaden 
the renormalized dispersion. We will obtain an analyti- 
cal formula for the spatiotemporal spectrum |afe(w)| 2 for 
the /3-FPU chain and compare the numerically measured 
width of the frequency peaks with the predicted width. 

In the Hamiltonian (|46[) . the nonlinear terms corre- 
sponding to the trivial resonances have been absorbed 
into the quadratic part via the effective renormalized dis- 
persion <Dfc. Therefore, the new effective Hamiltonian is 



N—l 

E 

k=l 



H = 



where 



ms ms 



N-l 

E 

k,l 1 m,s=l 



T ms ^ ms a%aia m a s , (55) 



3(3 



T. 



0. 



SNrj 2 
otherwise 



y/cjkUiu} m u> s , k 7^ m, and k ^ s 



(56) 



The new interaction coefficient Z*' a ensures that the 
terms that correspond to the interactions with trivial res- 
onances are not doubly counted in the Hamiltonian (|55|) 
and are removed from the quartic interaction. This new 
interactions in the quartic terms include the exact non- 
trivial resonant and non-trivial near-resonant as well as 
non-resonant interactions of the (2 — > 2)-type. 

We change the variables to the interaction picture by 
defining the corresponding variables b k via 

bk = ~a k e lGj *\ 

so that, the dynamics governed by the Hamiltonian (|55|) 
takes the familiar form 



ib k = 2 



N-l 

E 

Lm,s— 1 



rpkl Akl 1*7 I 11 



(57) 



where w^,. = Wfe+W;— ui TO — £> s [23j]. Without loss of gener- 
ality, we consider only the case of k < N/2. As we have 
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noted before, only for a very small number of quartets 
does vanish exactly, i.e., u)^ s = 0. We separate the 
terms on the RHS of Eq. (|57|) into two kinds — the first 
kind with u!™ s = that corresponds to exact non-trivial 
resonances, and the second kind that corresponds to non- 
trivial near-resonances and non-resonances. Since, in the 
summation, the first kind contains far fewer terms than 
the second kind, and all the terms are of the same or- 
der of magnitude, we will neglect the first kind in our 
analysis. Therefore, Eq. (|57|) becomes 



ib k 



JV-l 

E 



rpkl A kl 1*1 

1 ms LA ms°l °" 



(58) 



where the prime denotes the summation that neglects the 
exact non-trivial resonances. 

The problem of broadening of spectral peaks now be- 
comes the study of the frequency spectrum of the dy- 
namical variables bk(t) in thermal equilibrium. This is 
equivalent to study the two-point correlation in time of 
b k (t) 



C k (t) = (b k (t)bt(0)), 



(59) 



where the angular brackets denote the thermal average, 
since, by Wiener-Khinchin theorem, the frequency spec- 
trum 



(60) 



where J -1 is the inverse Fourier transform in time. Un- 
der the dynamics (155)) , time derivative of the two-point 
correlation becomes 



C k (t) 



= <-2i£ T^6r(t)6 ro (t)6 s (t)e-- t A« s 6U0)) 



l,m,s 



2i£ r£ 



kl „i<i£, s t jkl (4-\/\kl 
ms e J ms\ L ) L *ms' 



(61) 



where 



j^{t) = m)bi$)b m {t)b.(t))- 



In order to obtain a closed equation for Ck (t) , we need to 
study the evolution of the fourth order correlator J k „ (t) . 



We utilize the weak effective nonlinearity in Eq. (|55|) [121 ] 
as the small parameter in the following perturbation 
analysis and obtain a closure for Ck{t), similar to the 
traditional way of deriving kinetic equation, as in 0, l3l| . 
We note that the effective interactions of renormalized 
waves can be weak, as we have shown in [T3 |. even if 
the /3-FPU chain is in a strongly nonlinear regime. Our 
perturbation analysis is a multiple time-scale, statistical 
averaging method. Under the near-Gaussian assumption, 
which is applicable for the weakly nonlinear wave fields 
in thermal equilibrium, for the four-point correlator, we 
obtain 



j£(t)Al = C k (i)C,(0)(« + 8 k S l J. 



(62) 



Combining Eqs. (|56j) and (|62|) . we find that the right- 
hand side of Eq. (|5T)> vanishes because 



rpkl jkl /j.\\kl r> 



(63) 



Therefore, we need to proceed to the higher order con- 
tribution of Jms(t)- Taking its time derivative yields 



t kl <t)A 

J ms\ L )'-* 



kl 
ms 



([bt{t)b m (t)b s (t)+bt(t)b m (t)b a (t) 

&f(*)&m(*)6.(*K(0)>A«,. (64) 



Considering the right-hand side of Eq. (f64|) term by term, 
for the first term, we have 



(bt(t)b m (t)b s (t)bt(0))A 



kl 
ms 



[2i£ ffcb a (t)b}(t)b*(f) 



xe 



b m {t)b s (t)b%(0) )A 



kl 

ms " 



(65) 



We can use the near-Gaussian assumption to split the 
correlator of the sixth order in Eq. (f6"5|) into the product 
of three correlators of the second order, namely, 



(buo)b m ®b.(t)b a (t)bW)b;)K 

= Ck{t)n m n s 8 a (8™8^ + 5™5p 



kl 
m s 



where we have used that n r 
becomes 



C m (0). Then, Eq. 



(bt(t)b m (t)b s (t)b* k (0))A 
= Aif lk s C k (t)n m n s 



kl 
ins 



•A 



kl 

ms ' 



(66) 



Similarly, for the remaining two terms in Eq. (|64|) . we 
have 



kl 



{bf(t)b m (t)b s (t)b* k (0))A 
= -AifZt s Ck(t)nin s e lQ ™ A, 



kl 

ms > 



(67) 



and 



(b:(t)b m (t)b s (t)bl(0))A% s 
= -Aif^C k (t) ni n m e^ A^ 



(68) 

respectively. Combining Eqs. (|66|) . (f67|) . and (168)) with 
Eq. ([Ml), we obtain 



'^ms 

(t)A 

n 



kl 

m s 



ms 



Aif% s C k {t)e 
x (n m ii s - nin m - mn s ). (69) 



Equation (j6"9")) can be solved for J™ s it) under the assump- 
tion that the term e~ luJ ™ st oscillates much faster than 
Ck(t). We numerically verify [Fig. [5] below] the validity 
of this assumption of time-scale separation. Under this 
approximation, the solution of Eq. ([69j) becomes 



kl I 

ms 



■£',(*) A^, = 4T^C fc (t)A, 
x(n m n s - run 



1 



_ Qjkl 

nifis) 



(70) 
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Plugging Eq. ([70)) into Eq. (|6ip . we obtain the following 
equation for C k {t) 



\ 2 1 _ P iw„ s t 



C k {t) = 8iC k (t) V'( 

l,m,s 

x (n m n s - nin s - nin m ). 



(71) 



Since in the thermal equilibrium nu is known, i.e., n k = 
(\b k (t)\ 2 ) = e/Cj k [Eq. (25)], Eq. ((71]) becomes a closed 
equation for C/.(t). The solution of Eq. (JTTJ yields the 
autocorrelation function Ck(t) 



In 



C fc (0) 



i 



x(mn s 4- njn m - n rn n s )A 



kl 

ms " 



(72) 



Using this observation, together with Eq. (|56p. finally, we 
obtain for the thermalized /3-FPU chain 



In 



9fi 



2/12 



C k (t) _ 
C k (0) 8JVV 



■jjk 



ma 



\u)~,A i ,kl -i 



KL) 



(73) 



Equation (f?3")) gives a direct way of computing the cor- 
relation function of the renormalized waves a k , which, 
in turn, allows us to predict the spatiotcmporal spec- 
trum |S,fe(a;)| 2 . In Fig. [|3a), we plot the analytical pre- 
diction (via Eq. (|73p ) of the spatiotemporal spectrum 
\a k {uj)\ 2 = \b k {uj-Lo k )\ 2 = g-^C^^u-Uk). By compar- 
ing this plot with the one presented in Fig.[6](b), in which 
the corresponding numerically measured spatiotemporal 
spectrum is shown, it can be seen that the analytical 
prediction of the frequency spectrum via Eq. (|73p is in 
good qualitative agreement with the numerically mea- 
sured one. However, to obtain a more detailed compar- 
ison of the analytical prediction with the numerical ob- 
servation, we show, in Fig. [71 the numerical frequency 
spectra of selected wave modes with the corresponding 
analytical predictions. It can be clearly observed that 
the agreement is rather good. One of the important 
characteristics of the frequency spectrum is the width of 
the spectrum. We compute the width W(f) of the spec- 
trum f(uj) by 



W(f) 



max w f(u>) ' 



(74) 



In Fig. (8J we compare the width, as a function of the 
wave number k, of the frequency peaks from the nu- 
merical observation with that obtained from the analyt- 
ical predictions. We observe that, for weak nonlinearity 
{(3 = 0.125), the analytical prediction and the numeri- 
cal observation are in excellent agreement. In the weakly 
nonlinear regime, this agreement can be attributed to 
the validity of (i) the near-Gaussian assumption, and 




FIG. 6: (a) Plot of the analytical prediction for the spa- 
tiotemporal spectrum \&k(uj)\ 2 via Eq. ((73}. (b) Plot of the 
numerically measured spatiotemporal spectrum |afc(u)| 2 . The 
parameters in both plots were N = 256, (3 = 0.125, E = 100 
and T) = 1.06, 9 = 0.401. r\ and 6 were computed analyti- 
cally via Gibbs measure. The darker gray scale correspond to 
larger values of |afc(u;)| 2 in ui-k space. [max{ — 8, In |fifc(uj)| 2 } 
is plotted for clear presentation]. 
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FIG. 7: Temporal frequency spectrum |afc(w)| 2 for k — 30 (left 
peak) and k — 50 (right peak). The numerical spectrum is 
shown with pluses and the analytical prediction [via Eq. (|73|l ] 
is shown with solid line. The parameters were N = 256, 
(3 = 0.125, E = 100. 



(ii) the separation between the linear dispersion time 
scale and the time scale of the correlation C k (t). This 
separation was used in deriving the analytical predic- 
tion [Eq. (|73p ]. However, when the nonlinearity becomes 
larger (/? = 0.25 and (3 = 0.5), the discrepancy between 
the numerical measurements and the analytical predic- 
tion increases, as can be seen in Fig. [H Nevertheless, 
it is important to emphasize that, even for very strong 
nonlinearity, our prediction is still qualitatively valid, as 
seen in Fig. [5] In order to find out the effect of the umk- 
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FIG. 8: Frequency peak width W(\ak{u) | 2 ) as a function of 
the wave number fc. The analytical prediction via Eq. (|73(l 
is shown with a dashed line and the numerical observation is 
plotted with solid circles. The parameters were N = 256, E = 
100. The upper thick lines correspond to j3 — 0.5, the middle 
fine lines correspond to /3 = 0.25, and the lower solid circle 
and dashed line (almost overlap) correspond to /3 = 0.125. 



lapp scattering due to the finite size of the chain, we also 
computed the correlation [Eq. (|73[) ] with the "conven- 
tional" (^-function 5^ s (i.e., without taking into account 
the umklapp processes) instead of our "periodic" delta 
function A^ s . It turns out that the correlation time is 
approximately 30% larger if it is computed without umk- 
lapp processes taken into account for the case N = 256, 
j3 = 0.5, E = 100. It demonstrates that the influence 
of the non-trivial umklapp resonances is important and 
should be considered when one describes the dynamics of 
the finite length chain of particles. Finally, in Fig. [HI we 
verify the time scale separation assumption used in our 
derivation, i.e., the correlation time of the wave mode k 
is sufficiently larger than the corresponding linear disper- 
sion period tk = 2n /H>k- In the case of small nonlinearity 
(f3 = 0.125), the two-point correlation changes over much 
slower time scale than the corresponding linear oscilla- 
tions — the correlation time is nearly two orders of mag- 
nitude larger than the corresponding linear oscillations 
for weak nonlinearity f3 = 0.125, and nearly one order 
of magnitude larger than the corresponding linear oscil- 
lations for stronger nonlinearity (3 = 0.25 and (3 = 0.5. 
This demonstrates that the renormalized waves have long 
lifetimes, i.e., they are coherent over time-scales that are 
much longer than their oscillation time-scales. 



VII. CONCLUSIONS 

We have studied the statistical behavior of the nonlin- 
ear periodic lattice with the nearest neighbor interactions 
in thermal equilibrium. We have extended the notion of 
normal modes to the nonlinear system by showing that 
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FIG. 9: Ratio, as a function of k, of the correlation time Tk 
of the mode k to the corresponding linear period tk = 2n/il>k- 
Circles, squares, and diamonds represent the analytical pre- 
diction for /3 = 0.5, f3 — 0.25, and /3 = 0.125 respectively. 
Stars, pentagrams, and triangles correspond to the numerical 
observation for — 0.5, /3 — 0.25, and ft = 0.125 respectively. 
The parameters were N = 256, E = 100. The ratio is suffi- 
ciently large for all wave numbers k even for relatively large 
/3 = 0.5, which validates the time-scale separation assumption 
used in deriving Eq. (|70|l . The comparison also suggests that 
for smaller /3 the analytical prediction should be closer to the 
numerical observation, as is confirmed in Fig. [8] 



regardless of the strength of nonlinearity, the system in 
thermal equilibrium can still be effectively characterized 
by a complete set of renormalized waves, in the sense that 
those renormalized waves possess the Rayleigh- Jeans dis- 
tribution and vanishing correlations between different 
wave modes. In addition, we have studied the property 
of dispersion relation of the renormalized waves. The re- 
sults we obtained in Section [TT] are general and can be 
applied to the large class of nonlinear systems with the 
nearest neighbor interactions in thermal equilibrium. 

We have further focused our attention on the particu- 
lar system with the nearest neighbor interactions — the 
famous FPU chain. We have confirmed that the gen- 
eral renormalization framework that we discussed above 
is consistent with the numerical observations. In partic- 
ular, we have shown that the renormalized dispersion of 
the thermalized /3-FPU chain is in excellent agreement 
with the numerical one for a wide range of the nonlin- 
earity strength. We have further demonstrated that the 
renormalized dispersion is a direct consequence of the 
trivial resonant interactions of the renormalized waves. 
Using a self-consistency argument, we have found an ap- 
proximation of the renormalization factor via a mean- 
field approximation. In addition, we have used the mul- 
tiple time-scale, statistical averaging method to obtain 
the theoretical prediction of the spatiotemporal spectrum 
and demonstrated that the renormalized waves have long 
lifetimes. We note that the results obtained here can be 
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extended to general nonlinear potentials with the nearest 
neighbor interactions. 

The scenario of the wave behavior in the thermal 
equilibrium we obtained here may suggest a theoretical 
framework for the application of the wave turbulence to 
/3-FPU in the situation of near-equilibrium. 
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We start with the case of small nonlinearity /3 — > 0. Sup- 
pose in the first order of the small parameter (3 the tem- 
perature has the following form 

6(/3) = e +/39 u (A6) 

where O = O(l) and B\ = 0(1). We find the values of O 
and 6*i using the constraint (|A5|) . We use the following 
expansions in the small parameter (3 

e -^miy 2 +^) dy 

-co 

= yf^(4+(^-30 o )/3)+O(/? 2 ), (A7) 



APPENDIX A: LIMITING BEHAVIORS OF 77 
FOR THE THERM ALIZED /3-FPU CHAIN 



We change variables yj = qj 
nian for /3-FPU to obtain 



in the Hamilto- 



dy 



J^Vh (40o + (60i - 150g)/J) 



0(/3 2 ), (A8) 



H(p,y) 



N 

£ 

3=1 



2*' 



/? 4' 
4^ 



(Al) 



Next, we compute the pdf's for the momentum and dis- 
placement. Any pj is distributed with the Gaussian pdf 
z p = C p exp(— p 2 /2) and any yj is distributed with 
the pdf z y = C y exp(-0- 1 (y 2 + /3y 4 /2)/2), where C p and 
C v are the normalizing constants. As we have discussed, 
the renormalization factor r\ of the /3-FPU system in ther- 
mal equilibrium is given by Eq. (|25[) . and its approxi- 
mation via the self-consistency argument rj sc is given by 
Eq. ([52")) . Here, we compare the behavior of both formulas 
in two limiting cases, i.e., the case of small nonlinearity 
(3 — > and the case of strong nonlinearity (3 — > 00. We 
will use the following expressions for the average density 
of kinetic, quadratic potential and quartic potential parts 
of the total energy of the system 



(K) 


1 ^ 


N 


N ' 




3 


(U) 


1 ^ 


N 


N 




3 


(V) 


(3 


N 


4iV 



3=1 



(Vl 



1 



(y 2 



N 

3=1 



(A2) 
(A3) 
(A4) 



In a canonical ensemble, the temperature of a system is 
given by the temperature of the heat bath. By identifying 
the average energy density of the system with e = E/N 
in our simulation (a microcanonical ensemble), we can 
determine as a function of e and (3 by the following 
equation 



1 

N 



((K) + (U) + (V)) 



e. 



(A5) 



{y 2 + | 2/ 4 )e -^ 3 y(^+^) dy 



% (40 o + (60! - 90 2 )/3) + 0(/3 2 ). (A9) 

Then, in the first order in /3, Eq. (|A5|) becomes 
/ |v / 0o"(40 o + (60i- 90 2 )/3) 



/30! 



fV^(4+(|f-30o)/3 



= 2e, 



and we obtain 8q = e and 0i = (3/4)e 2 . Therefore, for 
the average kinetic energy density, we have 



N 2* 



(Aio) 



and, for the average quadratic potential energy density, 
we have 



(CO 1 y/fv^ (40q + (60i - 150 2 )/3) 
N ~ ~ 2 yfV^(4+( 2 e-30 o )/3) 

= ^-?e 2 /3 + 0(/3 2 ). 
Finally, we obtain Eq. (fS"3"]) . i.e., for small /3 



ry=l + -e/3 + 0(/3 2 ). 



(AH) 



(A12) 



Similarly, from Eq. (|52j) . we find the small /3 limit of the 
approximation rj sc 

Vs C = l + hl3 + 0((3 2 ). 

Now, we consider the case of strong nonlinearity f3 — > 
00. From Eq. (| A5|) . we conclude that temperature in the 
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large (3 limit, which we denote as 9oo, stays bounded, i.e., 
< 6*oo < 2e, and, in the limit of large (3, we obtain for 
Eq. JAU) 



J: 



r y4 dy 



2e. 



(A13) 



After performing the integration, we obtain 9^ = (4/3)e, 
and the average kinetic energy density becomes (K) /N = 
(2/3)e. For the average quadratic potential energy den- 
sity, we have 



(U)_ _ l QVg dy = £(|) /4e 
2 /T^e-^dy r (i) V3/3 



V 



^ ( A14 ) 



For the renormalization factor, we obtain the following 
large (3 scaling 



Similarly, for the approximation of rj sc , we obtain A — 
Cy/&P, B = 4ej3, and C = 2v / 3T(3/4)/r(l/4). There- 
fore, the large (3 scaling of rj sc becomes 



C + VC 2 + 16 
Vsc = \ 7; e*/3* 



(A16) 



r(l) 



Var(i) 



— e 4 /3< 



(A15) which yields Eq. ([51)1 . 
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